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ABSTRACT 

We investigate the effects of the change of cosmological parameters and star formation 
(SF) models on the cosmic SF history using cosmological smoothed particle hydrody- 
namics (SPH) simulations based on the cold dark matter (CDM) model. We vary the 
cosmological parameters within 1-cr error from the WMAP best-fit parameters, and 
find that such changes in cosmological parameters mostly affect the amplitude of the 
cosmic SF history. At high redshift (hereafter high-z), the star formation rate (SFR) is 
sensitive to the cosmological parameters that control the small-scale power of the pri- 
mordial power spectrum, while the cosmic matter content becomes important at lower 
rcdshifts. We also test two new SF models: 1) the 'Pressure' model based on the work 
by Schaye & Dalla Vecchia (2008), and 2) the 'Blitz' model that takes the effect of 
molecular hydrogen formation into account, based on the work by Blitz & Rosolowsky 
(2006). Compared to the previous conventional SF model, the Pressure model reduces 
the SFR in low-density regions and shows better agreement with the observations of 
the Kennicutt-Schmidt law. This model also suppresses the early star formation and 
shifts the peak of the cosmic SF history toward lower redshift, more consistently with 
the recent observational estimates of cosmic SFR density. The simulations with the 
new SF model also predict lower global stellar mass densities at high-z, larger popu- 
lations of low-mass galaxies and a higher gas fraction in high-z galaxies. Our results 
suggest that there is room left in the model uncertainties to reconcile the discrepancy 
that was found between the theory and observations of cosmic SF history and stellar 
mass density. Nevertheless, our simulations still predict higher stellar mass densities 
than most of the observational estimates. 

Key words: method : numerical — galaxies : evolution — galaxies : formation — 
galaxies : high redshift — galaxies : mass function — cosmology : theory 



1 INTRODUCTION 

Star formation is a fundamental and important physical 
process for understanding galaxy formation. Stars emit a 
large fraction of the observed photons, and their formation 
and evolution significantly influence galaxy evolution. Over 
the past decade, various high-2 observations started to 
unveil the rough h i story of cosmic s tar formation (e.g. , 
Madau et all Il996l: iLillv et al.l Il996l: ISteidel et ai] Il996l: 



Dickinson et all l2003l: iGiavalisco et all l2o61" lOuchi et all 
2004 iHopkins fc BeacomI l2006l : iBouwens et~ai] l2007l : 
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lYan et all l2008l : iBouwens etHI |2009| . and references 
therein). There are yet considerable uncertainties in the 
dust extinction correction, the faint-end slope of the 
luminosity function, and the stellar initial mass function 
(IMF). However, these observations show a rough picture 
of the cosmic SFR density (p*), which gradually increases 
from high-z to z ~ 4, peaks at z ~ 2 — 4, and then rapidly 
declines from z ~ 2 to z = 0. It is then of significant interest 
to see if the standard model of cosmic structure formation 
based on the CDM model can reproduce the observed 
features of cosmic SF history, and to further develop insight 
into the physical processes which govern the cosmic star 
formation. 

Numerical simulations have been use d to investigate 
the c osmic star forma t ion history (e.g., ICen fc Ostriker 
1 19921 : iKatz et al.lll996al : lYepes et all 1 19971 ; iNagamine et al 
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2OOOI; Weinberg et al. 


2002 


: SDringel & Hernauistl l2003bl: 


NaEamine et alj 2004,120061: 


Davell2008l), however, most the- 



oretical predictions remain somewhat uncertain due to the 
following reasons. First of all, the physics of star formation 
and its feedback have not been clearly understood yet, be- 
cause they are inherently complex problems, involving non- 
linear dynamics, ra diative and chemical pro cesses on wide 
ranee of scales fe.g.. lMcKee fc Ostrikeij|2007l ). Secondly, the 
resolution limitation of simulations forbids the detailed mod- 
eling of SF and its feedback from first principles, and cos- 
mological simulations must rely on empirical SF laws (e.g, 
ISchmidt|[l959l : lKennicutt|[l998al lbh. The simulation results 
could depend on the details of the adopted models of SF 
and its feedback, and there are still significant freedoms in 
the formulation of these models. 

Originally, numerical cosmologists studied the cosmic 
SF history with the motivation to use it as a probe of cos- 
mological models, and examine the effects of largely differ- 
ent cosmological mod els on (e.g., iNagamine et al.|[200ol : 
[Weinberg et alll2002l ). This is because different cosmologi- 
cal models result in different power spectra of matter density 
fluctuations, which would be reflected as different shapes of 

Since then, the situation has dramatically changed, and 
we now know the val ues of our cos mological |)arameters 
within ~10% accuracy (jKomatsu et a l. 2009, 2010), and the 
resolution of numerical simulations has also improved signif- 
icantly since late 1990's. Back then, cosmological hydro sim- 
ulations with 64^ particles were standard, but now we can 
go up to 300'^ — 400^ particles. Given these advancements, 
it would be worthwhile to quantify the changes in in 
more detail when we deviate the cosmologica l parameters by 
1-cr error f rom the WMAP best-fit values (jKomatsu et al.l 
|2009| , |2010| ). An interesting recent change is that the values 
of (78 and the power-index of primordial power spectrum Us 
have become smaller with as < 0.9 and Us < 1.0. The goal 
of this paper is to quantify the dependence of Pi,{z) on cos- 
mological parameters and the adopted SF model in detail, 
as it has not been documented in the literature yet. 

A very comprehensive work was performed by 
ISpringel fc Hernauistl (|2003bl , hereafter SH), in which they 
estimated ^^(z) using a large set of cosmological SPH simu- 
lations with a n ovel subgrid model for multiph ase interstellar 
medium (ISM) (|Springel fc Hernauistll2003al '). The parame- 
ters in their SF model were con strained by the observed 
Kennicutt- Schmidt law (jKennicutt 1998ai)bi). They clearly 
showed how p*(z) depends on the numerical resolution, but 
they did not examine the eflects of different cosmologies and 
SF models in detail. 

Based on the work of SH, iHernauist fc Springell l|2003l , 
hereafter HS model) provided a simple, analytic reasoning to 
identify physical processes that drive the evolution of p-t(z). 
They argued that the early phase of cosmic star formation 
is driven by the gravitational structure formation, and pro- 
vided an analytic formula for p*(z), which includes the Hub- 
ble parameter H{z) and the SFR function as a function of 
halo mass. Compared to the observational estimates, the HS 
model with the WMAP cosmology predicts a higher at 
high-z with a peak at z « 5 — 6, resulting in higher stellar 
mass densities ( i i*) at z > 1 than the observations. Later, 
INagamine et al.l (|2004l ) compared two different kinds of nu- 
merical simulations (Eulerian TVD and SPH) with obser- 



vations, and found that the numerical simulations generally 
predicted higher SFR at high- 2:. 

Although the SH model provides a novel description 
of star formation and its feedback, this model is still not 
perfect. As mentioned above, the star formation in the 
SH model is constrained by the Kennicutt-Schmidt law 
which is the correlation between the average star forma- 
tion surface density (Esfr) and average total gas surface 
density (Egas). Recent observations, which spatially resolve 
star formation in a given galaxy, show that EsFH,(r) is 
a function of the EH2(r), and Egas(r)/EH2 (r) is not con- 
stant. This suggests that Esfr is indeed a function of 
molecular hydrogen surface density (Eh,), rather than E| 



(Wong fc Blitz 2002: He ver et al.ll2004l : iBlitz fc Rosolowskvl 
^2006i : iBigiel et al. 2008). Therefore, it will be important to 
take into account the effect of H2 density in our SF model. 
In this paper, we examine new SF models which consider 
the contribution of H2, and test the effects on p*(z) and 
star formation on small scales. 

This paper is organized as follows. In Section [2l we de- 
scribe our simulations with a focus on the new SF models. 
In Section [S] we present the effects of differing cosmologi- 
cal parameters on the cosmic SF history. In Section [4] we 
compare the two new SF models with the SH model. The 
effect of new SF models on Pic{z) and galaxy formation is 
investigated in Section [S] Finally we summarize our flndings 
in Section |6l 



2 NUMERICAL TECHNIQUE 

We use the updated version of the Tree-particle-mesh 
(TreePM) s moothed parti cle hydrodynamics (SPH) code 
GADGET-2 l|Springell2005l ). Our conventional code includes 
radiative cooling by H, He, and metals (jChoi fc Nagamin^ 
200!t) , heating by a uniform UV back ground of a modi- 
ficd Haar dt fc Madaul (| 19961) spectrum ijKatz et all Il996bl : 
iDavc et al.l Il999l l. star formation, supernova feedback, a 
phenomenological model for galacti c winds, and a sub- 
resolution model of multiphase ISM jSpringel fc Hernauistl 
2003a). In this multiphase ISM model, the high-density ISM 
is pictured to be a two-phase fluid consisting of cold clouds in 
pressure equilibrium with a hot ambient phase. Cold clouds 
grow by radiative cooling out of the hot medium, and this 
material forms the reservoir of baryons available for star for- 
mation. Since the details of the t r eatm ent are described by 
both SH and I Choi fc Nagamind (|2009l ). we do not repeat 
them here. 



2.1 Cosmological parameters 

In this paper, we adopt the following flducial cosmology 
which is consistent with the latest WMAP result: Qm = 
0.26, = 0.74, Qb = 0.044, h = 0.72, = 0.96, and 
(Tg — 0.80. In order to see the effect of each cosmologi- 
cal parameter on pt,{z), we vary one parameter at a time 
from the flducial model while the other parameters are kept 
flxed. Except erg, we vary each parameter from the fiducial 
value by approximately one sigma error in the WMAP result 
(|Komatsu et al.ll2009l . I2OI0I ). The variation of ag is chosen 
to be 0.1 to include the previously popular ag = 0.9 in our 
parameter set. The tested cosmological models are in the 
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range of Us = 0.92 - 1.0, as = 0.7 - 0.9 ^6 = 0.038 - 0.05, 
{Qm, ^Ia) = (0.24, 0.76) and (0.28, 0.72) Q. The runs used in 
this comparison are hsted in Table [1] 

Star formation at high redshift takes place mostly in 
low-mass galaxies, and it progressively shifts to massive sys- 
tems at lower redshifts. We need a high-resolution simula- 
tion to resolve the star formation in low-mass galaxies and 
a large box size to include a large number of massive galax- 
ies, which requires substantial computational resources. In 
order to get around this difficulty, we employ a large num- 
ber of runs with different resolution and volumes. Because 
of this, we limit the particle number of our fiducial runs to 
2 X 216'' (gas -I- dark matter) particles in a comoving box of 
(10/i"^Mpc)^ (hereafter the N216L10 series). We stop this 
series at z — 2.75, as it misses the long wavelength pertur- 
bations at lower redshifts. The resolution of N216L10 series 
is not adequate to properly simulate the entire cosmic SF 
history, but it is sufficient to show the differences due to the 
variation of the cosmological parameters. 



2.2 Star formation models 
• SH model: 

In the multiphase ISM model of SH, the star formation is 
modelled as follows. If the dense gas is Jeans unstable and is 
rapidly cooling, a fraction of the gas mass is converted into 
a star particle with a rate 



p. = (1 - P)pc/ts 



(1) 



where pc is the gas density of the cold cloud, and /3 is the 
mass fraction of high-mass stars that instantly die as su- 
pernovae, determined by the stellar initial mass function. 
The star formation time-scale tsFR is taken to be propor- 
tional to the local dynamical time of the gas: tsFnip) = 
to {p/ pth)~^^^ , where the value of = 2.1 Gyr is chosen 
in isolated disk galaxy simulations to match the Kennicutt- 
Schmidt law: 



SsFR ~ 





A(Egas/lM0 pc- 



if Eg 
if Eg 



< Eth 
> Eth, 



(2) 



where Eth is the SF threshold surface density. Observations 
suggest that A = 2.5 ± 0.7 yr~^ kpc~ ^, n = 1.4 ± 0.15, 
and Eth ~ lOMopc"^ l|Kennicuttlll998al lbf). 

The supernova explosions add thermal energy to the 
hot phase of the ISM and evaporate cold clouds. This 
is described as pc = Cppc/tt,, where the feedback effi- 
ciency parameter "C" has the density dependence C{p) = 
Co (p/pth) following iMcKee fc OstrikeJ (|l977l '). This 
evaporation process of cold clouds establishes a tight self- 
regulation mechanism for star formation in the ISM, where 
the ambient hot medium quickly evolves toward an equilib- 
rium temperature. 

Although the SH model has shown significant improve- 
ments over the earlier SF models, it is still incomplete. We 
know from observations that star formation takes place in 
the molecular clouds. Therefore the next natural step is to 
implement alternative SF models which include the effect 
of molecular hydrogen (H2), and study the consequences on 



^ Tlie flat universe is one of the most robust constraint of WMAP, 
tlierefore we fix Clrn + f^A = 1 for this variation. 



the cosmic star formation history (e.g.. iGnedin et al]|2009l V 
Unfortunately, current cosmological simulations cannot re- 
solve the detailed structure of molecular clouds while solving 
the formation of thousands of galaxies on a >10Mpc scale, 
so we need a model to estimate the H2 fraction from the 
total gas density. Here, we consider two new SF models and 
compare them with the original SH model. 

• Blitz model: 

This mo del is based on the H 2 dens ity-pressure relation de- 
rived by iBlitz fc Rosolowskvl (|2006l ). They argued that the 
mean ratio of molecular to atomic hydrogen surface den- 
sity is related to the interstellar gas pressure as follows: 
Ehs/Ehi = (Pext/^li)", where Pext is the interstellar gas 
pressure, Po/k = 4.3 x 10* cm~^ K, and a x 0.92. Using this 
relationship, we can compute the amount of H2 from the 
total gas density and pressure. The projected SFR density 
in this model is 



EspR = O.leEg 



1 + 



Paxt 



(3) 



where e « 10 - 13 Gyr" ^ . iBlitz fc Rosolowskvl (|2006l ) argued 
that this relationshi p recovers the obse rved SFR surface den- 
sity better than the lKennicut j l|l998ah law, especially in the 
low density regime, where the molecular density is lower 
than the Hi density and Pext < Po- 

We assume that we can replace the surface density with 
the 3-dimensional density, and rewrite the above equation as 



1 + 



Pext 

"pT 



Gyr- 



(4) 



after adopting e = 10 Gyr ^ as suggested in 
iBlitz fc Rosolowskvl (|2006l l. We note that this SF law 



is very similar to the one adopted by iKravtsqy (|2003|) 
excep t for the pressure term in the bracket. iKravtsovl 
(j2003h assumed p* = pgas/r* with r* = 4 Gyr. 

When Pext > Po, the molecular surface density be- 
comes greater than the atomic hydrogen surface density. In 
this high-pressure regime, we assume that the SF law re- 
verts to the Kennicutt-Schmidt law, and adopt equation ((5}. 
We apply equation ||4]) only in the low-pressure regime with 
Pext < Po- We call this new model the 'Blitz' model. 

• Pressure model: 

This model explicitly formulates the conversion between gas 
surface density (Egas) and gas volume density (pgas). Pre- 
viously, we assumed Egas/EsFR = Pgas/p*, which is only 
true if the disk scale-height is consta nt or t h e equ ation 
of state (EOS) behave s as P oc p^. ISchavd (|200ll ) and 
ISchave fc Dalla Vecchial l|2008l ) proposed the "Jeans column 
density" , and argued that the scale-height will be of the or- 
der of local Jeans scale for self-gravitating discs, because the 
density typically fluctuates on the local Jeans scale. 

One may argue that this model takes into account the 
effect of H2 cooling better than the previous SF models 
based on the three-dimensional gas density, because the disk 
instability leads to the collapse of molecular clouds and star 
formation occurs within them due to the H2 cooling. Ow- 
ing to the current limitation in computational power, it is 
still impossible to simulate the formation of molecular clouds 
explicitly in cosmological simulations. Therefore considering 
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Runs 








ris 


0-8 


"IDM 




FiQucial 


v. ZD 


n 7/1 

U. /4 




u.yo 


U.o 


5.96 


X 


10^ 


1.21 


X 


lO'' 


High ris 


0.26 


0.74 


0.044 


1.00 


0.8 


5.96 


X 


10** 


1.21 


X 


10« 


Low rig 


0.26 


0.74 


0.044 


0.92 


0.8 


5.96 


X 


10^ 


1.21 


X 


10^ 


High as 


0.26 


0.74 


0.044 


0.96 


0.9 


5.96 


X 


10^ 


1.21 


X 


10" 


Low <T8 


0.26 


0.74 


0.044 


0.96 


0.7 


5.96 


X 


lO'' 


1.21 


X 


lO'' 


High Qrn 


0.24 


0.76 


0.044 


0.96 


0.8 


5.41 


X 


10^ 


1.21 


X 


lO'' 


Low Qtji 


0.28 


0.72 


0.044 


0.96 


0.8 


6.51 


X 


10^ 


1.21 


X 


10« 


High 


0.26 


0.74 


0.038 


0.96 


0.8 


6.12 


X 


10'^ 


1.05 


X 


10" 


Low Qf, 


0.26 


0.74 


0.050 


0.96 


0.8 


5.79 


X 


10*^ 


1.38 


X 


10" 



Table 1. The simulations employed in the cosmological parameter test. All the simulations are N216L10 series whose general properties 
are shown in Table [2] The mass resolution varies slightly due to the variation of the cosmic matter content 



Series 


Box-size 






m,DM 


m 


-gas 


€ 


■^cnd 


N216L10 


10.0 


2 


X 216^ 


5.96 X 10" 


1.21 


X 


10" 


1.85 


2.75 


N400L34 


33.75 


2 


X 400^ 


3.49 X lO'^ 


7.31 


X 


10" 


3.375 


1.0 


N400L100 


100.0 


2 


X 400^ 


9.12 X 10* 


1.91 


X 


10* 


6.45 


0.0 



Table 2. The two series of simulations employed for the comparison study of star formation models. The box-size is given in units of 
/i^^Mpc, A^p is the particle number of dark matter and initial gas (hence X 2), mj^M and rrigas are the masses of dark matter and gas 
particles in units of h~^MQ, respectively, e is the comoving gravitational softening length in units of /i~^kpc, and Zg^d is the ending 
redshift of the simulation. The value of e is a measure of spatial resolution. Note that the star particle mass is a half of the initial gas 
particle mass. We implement three star formation models (the SH model, the Blitz model, and the Pressure model) for the N216L10 
series and two star formation models (the SH model and the Pressure model) for the N400L100 series. 



the Jeans instability in a disk could be a useful approxima- 
tion to take account of the effect of disk instability, which 
leads to the molecular cloud formation, in cosmological sim- 
ulations. 

Using the Jeans column density surface and vol- 

ume gas densities are related as follows: 



Egas ~ Sg,.J = PgasLj — ^ —fgPtot, (5) 

where Lj = Caj \fGptot is the Jeans length, Ca = 
\fyPtotlp^s is the local sound speed, is the mass frac- 
tion in gas (i.e., pgas = fgPtot), and fth is the fraction of 
mid-plane pressure that is thermal (i.e., P = fthPtot, where 
P is the thermal pressure). 

Throughout our calculation, we assume /<, — fth, there- 
fore fgPtot = {fg/fth)P = P- From equations ^ and ©, 
we can derive the new SF time-scale: 

W = |^=A-(lMePc-)"ap)'^-"'^\ (6) 

where we adopt n = 1.4 and 7 = 5/3 as the default values. 
In this SF model, we keep the multiphase ISM model of SH. 

We compute the SF threshold density as pth = 
G / fg{T,th/ ■ To compute Cs for the gas at the thresh- 
old density, we need to know the temperature and the mean 
atomic weight of the gas. Here, we assume Tth = SOOif and 
Pth = 1.2, because star formation occurs in the cold gas. 
The resulting number density is nth ~ 1.2h? cvcC^ , which 
is a factor of three higher than the original value used in 
the SH model. As we will discuss later, our new value of 
nth reproduces the proper threshold column density of star 
formation, as observed bv lKennicut"^ (|l998al ). 



The shape of the EOS is closely related to the feedback 
model. In the multiphase prescription of SH, the thermal 
energy from SN feedback pressurizes the ISM and makes 
the EOS steeper. iRobertson et al. (2004) later derived a fit- 
ting formula for this steep, effective EOS of the SH model. 
Although it is possible to use the effective EOS in the sim- 
ulation for the star-forming regions, we prefer to adopt a 
polytropic EOS for the star-forming gas to prevent an arti- 
ficial fragmentation discussed below. Note that we still apply 
the orig inal EOS from the SH model fo r the non-star form- 
ing gas. ISchave k, Dalla Vecchial l|2008l ) discussed the use of 
a polytropic EOS, 

P = KpllV. (7) 

Note that the effective polytropic index 7e/ / is not the same 
as the usual adiabatic index 7. If 7^/ / =4/3, the Jeans mass 
is independent of the gas density. If the Jeans mass decreases 
with density, it will lea d to the artificial fragmentation of gas 
(|Bate fc Burkertlll997l '). Therefore a polytropic index 7e// = 
4/3 prevents the artificial fragmentation while allowing the 
collapse to proceed, and we adopt this value as our default. 
The value of K in equation (O is computed by inserting this 
equation into equation ((S)), K = (G/7)pj^^°^^ Ej^j. 

With the new SF time-scale, new density threshold, and 
new EOS for the star-forming gas, we go back to equation ([l]) 
and complete the new SF model. We call this model the 
'Pressure' model. In Sectional we will compare pi,{z) from 
these three SF models using a series of simulation sets pre- 
sented in Table [2I 
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Figure 1. The cosmic star formation iiistory from tiie N216L10 series of simulations with different cosmological parameters (see Table[T|. 
The shape of the cosmic star formation history is not changed by the different cosmological parameters, but the amplitude of the star 
formation history is significantly affected. 



3 VARYING THE COSMOLOGICAL 
PARAMETERS 

Figure [T] shows the cosmic star formation histories of the 
N216L10 series with different cosmological parameters. The 
general shape of p*(2) is preserved when the cosmological 
parameters are varied, but the amplitude changes. Figure [1] 
shows that the effect of the primordial power spectrum on 
Pi is more significant than that of the global mass content of 
the Universe. The variance of Pi, is quite large at high red- 
shift when Us and ug are varied. This is because the values of 
Us and (Jg determine the primordial power spectrum, which 
governs the early structure formation. The global matter 
content (fi) also changes the power spectrum, but its con- 
tribution is mostly on the transfer function and the growth 
factor, whose influence becomes more important at lower 



redshifts. Therefore, it is expected that the effect of change 
in Tis and erg is more significant than that of Q. at high-z. 

Figure [1] also shows that the variance of p* becomes 
smaller at lower redshifts when ris and o-g are varied. At 
early times, the cosmic density is high, and the cooling time 
is short, therefore most dark matter halos can cool the gas 
and form stars. During this epoch, the cosmic SFR is mostly 
driven by the gravitational growth of dark matter halos. 
At lower redshifts, the cooling time becomes longer due to 
decreasing cosmic density and hotter intergalactic medium 
(IGM). Then the effect of gas physics starts to play a more 
important role in determining the SFR, and the tight cou- 
pling between dark matter halo growth and SFR becomes 
weaker. 

The cosmological matter densities (f^A, f^m, and Q.h) 
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Figure 2. Galaxy stellar mass functions at z = 3 for the runs with different cosmological parameters (see Table [TJ. The shape of the 
mass functions is hardly affected by the variation of cosmological parameters, while the amplitude of them is slightly changed by the 
variations. 



contribute to the cosmic SF history in two difTerent ways. 
The values of Ha and Qrn control the redshift evolution of 
the Hubble parameter, H{z), which affects the linear growth 
factor of structure formation and the expansion rate of the 
universe. Therefore the changes in SIa and Qm affect the 
p*{z) through the growth factor. At low redshift, the de- 
crease of the mean density of the Universe results in th e 
declining efficiency of gas cooling (|White fc FrenkI llQQlT ) . 
Hence, JIa and Qm continuously influence the cosmic star 
formation history, which is confirmed in Figure [T] 

The variation of Qi, simply changes the gas density for 
a given dark matter halo. The increase (decrease) of the 
gas density enhances (reduces) gas cooling. Therefore the 
effect of Qb becomes important at lower redshifts when the 
effect of gas cooling time is more important. Our simulations 



confirm that the cosmic SF history is mostly determined 
by the structure formation at high redshi ft and by the gas 
cooling at low redshift, as discussed bv iHernguist fc Springell 
(|2003h . 

Figure [2] shows the galaxy stellar mass functions at 
z — 3 for the runs with different cosmological parameters. 
The galaxies in the simulation were ide ntified using a simpli- 
fied variant of the SUBF IND algorithm (|Springel et al.ll200ll : 
IChoi fc Nagaminell2009l ). In this paper, we set the mass limit 
of the simulation galax;ies 32 stellar particle mass (see Ta- 
bleful. The shape of mass functions is hardly affected by the 
variation of the cosmological parameters, while the ampli- 
tude of them is slightly changed. One noticeable feature in 
Figure [2] comparing with Figure[T]is that the change of mass 
function owing to the variation of Us is relatively small. At 
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Figure 3. The projected SFR as a function of H I column density at ^ = 3 for the N216L10 scries with different cosmological parameters. 
The red line shows the prediction of the Kennicutt-Schmidt law. The cosmological parameters hardly change the shape of this relation. 
The blue dashed line is the current observational limit of the low SFR at 2: = 0. 



first glance, it might appear to be conflicting with the sig- 
nificant change in Pt{z). However, the SFR in Figure [T] is 
in a logarithmic scale, therefore the seemingly large relative 
difference in pi,{z) at high-z only results in a small difference 
in the normalization of the galaxy stellar mass function. In 
addition, the difference in p+(z) caused by ris becomes neg- 
ligible at z=3. Generally the qualitative trend in Figures [T] 
and [5] is the same, with the green dashed line at the top, 
the red curve in the middle, and the blue dotted line at the 
bottom. 

The cosmological parameters do affect the cosmic star 
formation history, but they hardly change the SF efficiency 
as a function of gas column density. Figure [3] compares 
the SFR in our simulations with the empirical Kennicutt- 
Schmidt law (red solid lines, with upper and lower lines 



showing the errors in Equation ([2])). All simulations show 
similar results and relatively good agreement with the 
Kennicutt-Schmidt law (except for the low column density 
end, as we will discuss later). Figure [3] represents the local 
star formation relation, while Figure [T] represents the global 
star formation. This comparison shows that, although the 
cosmological parameters infiuence the amplitude of cosmic 
SFR, they hardly alter the local SFR distribution. The lo- 
cal star formation is mostly controlled by the adopted SF 
model. In the next section, we will compare the effect of 
different SF models. 
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Figure 4. The projected SFR as a function of Hi column density at z = 3 for the N216L10 series with different SF models. The red 
lines show the empirical Kennicutt-Schmidt law in Equation l(2]l. All three models show reasonably good agreement with the Kennicutt- 
Schmidt law, however there are notable differences in the three results as we discuss in the text. The blue dashed line is the current 
observational limit of the low SFR at z = 0. The yellow shading in e ach panel indicates the range of observed data for 2 = galaxies 
taken from the locus of data points in Figure 15 of iBigiel et al] j200^ . 



4 COMPARING DIFFERENT STAR 
FORMATION MODELS 

In Section 12.21 we proposed two new SF models: the Blitz 
model and the Pressure model. In this section, we compare 
the results of these models with those of the original SH 
model. Figure |4] shows the projected SFR as a function of 
Hi column density, whic h we computed followin g the pre- 
scriptions as described in iNagamine et al] (120041 ) . The size 
of the projected pixels is equal to the gravitational soften- 
ing length in the simulation. The three SF models generally 
agree with the Kennicutt-Schmidt law, but there are some 
notable differences between the models particularly in the 
low density and low Ssfr regions, which we now discuss. 



The projected SFR from simulation is for 2 = 3 result, 
while the observed data is for z = 0. Most high-2; galax- 
ies are currently not resolved well, therefore it is difficult 
to measure the Kennicutt-Schmidt law for high-z galaxies. 
There are some ind irect inferences, such as the work by 
I Wolfe fc Ch^l|2006D . that suggest lower star formation effi- 
ciencies at high-redshift, but the uncer tainties are still very 
large. On the contrary, another work bv lBouche et al.l (|2007l ) 
actually suggests the opposite with a four times more effi- 
cient star-formation for z ~ 2 star-forming galaxies. There- 
fore, at this moment, it would be best to compare our sim- 
ulation results against the most robust local observations of 
the Kennicutt-Schmidt law. 
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First, all simulations including those in Figure [3] show 
a large population of low Sspr with log Esfr < —5 and 
log A'hi ~ 18 — 19. These low Esfr SF regio ns are below 
the lo wer limit of the current observations (e.g.. lBigiel et all 
|2008| ). which is indicated by the dashed horizontal line. The 
fraction of star formation occurring in these low Esfr re- 
gions are 0.017, 2.9, and 0.0052% for the SH, Blitz, and 
Pressure model, respectively, while the number fractions of 
columns in these low Esfr regions are 18, 72, and 31% for 
the SH, Blitz, and Pressure model, respectively. The total 
SFR in these low Esfr regions is negligible, because the 
largest contribution comes from high density and high SFR 
regions. However, there is a significant number of columns 
with these low Esfr. It is possible that the current observa- 
tions are missing the above fractions of SFR and the number 
of star-forming columns in the Universe. These low Esfr re- 
gions might be related to the outskirts of the disk, or the 
destruction of dwarf galaxies and the formation of tidal tails. 
It would be an interesting future topic of research to further 
investigate these low Esfr regions in our simulations. 

Second, the SH model overpredicts the Esfr at 
logA'ni J; 20.5, which is also true for all the runs shown 
in Fig ure El This feature was also noted byll agamine et al.l 
l|2004 ). This overprediction of Esfr is absent in the Blitz 
model and the Pressure model. The Blitz model considers 
the effect of H2 on the SFR in the low-pressure regions with 
low SFRs, and the lower SFRs in the low A'^hi regions are ex- 
pected in this model. This appea rs to be a favorable improve - 
ment of the SF model, because iBlitz fc Rosolowskvl (|2006l ) 
compared the observed Esfr against their SF model and 
the Kennicutt-Schmidt law for several local galaxies, and 
found that their model shows better agreement with the ob- 
servations than the Kennicutt-Schmidt law at the outskirts 
of disks, where the gas density is lower. In our cosmolog- 
ical simulations, the Blitz model gives overall lower Esfr 
than the observed Kennicutt-Schmidt law (but still within 
the observed range shown by the yellow shading) as shown 
in the top right panel of Figure U There might still be some 
room for improvement of this model in our simulations. 

Third, the Pressure model shows more favorable fea- 
tures compared to both the SH and the Blitz model (see 
Figure [4}. As described in Section [2.21 the Pressure model 
adopts a higher SF threshold density than the SH model, 
and this has two effects. One is the reduction of Esfr in 
low A^Hi regions compared to the SH model, and the other 
is the shift of the SF cut-off column density towards higher 
value at around log A^hi ~ 19.5. Since the current observa- 
tions only probe down to log Esfr — —4.5, the turn-down 
may appear to occur at log Esfr ~ 20 — 21 fe.g.. lKennicuttl 
ll998al : lBigier et al.|[2008l ). The changes seen in the Pressure 
model are favorable improvements over the SH and Blitz 
models, as they bring the simulation results closer to the 
observations. In addition, the Esfr for a given A'^ei is lower 
than the Kennicutt-Schmidt law in the Blitz model, which 
makes the Blitz model somewhat unfavorable. This may be 
solved by increasing the time-scale parameter e in the model, 
however, here we choose not to change the parameter values 
suggested by the original authors, because they determined 
those values by observations. For these reasons, we favor the 
Pressure model over the SH and Blitz models. 

Lastly, the slope of Esfr in the Pressure model is 
slightly shallower than the SH and the Blitz model at high- 




SH model 
Pressure model 
HS model 



Figure 5. The cosmic star formation history from our simu- 
lations with different SF models. The solid lines are from the 
N216L10 series, and the dashed lines are from the N400L100 
series. The N216L10 series represent the high-z SFR better, 
and the N400L100 series represent the low-z SFR better. We 
compare our results with the previous theoretical model of 
iHernouist &: Springell ||2003| . the HS model; blue long-dashed 
line). The cy a n sha ding is the observed range of SF history from 
iKistler et al. I l l200gD . Th e yellow shading is the locus of the ob- 
served data compiled bv lNagamine et al.l ||2006| V Both compila- 
tions of data considered the dust extinction correction. This figure 
shows that the peak of the SFR density shifts to a lower redshift 
in the Pressure model compared to the SH model. 



A'hi with n ~ 1.3 (see Equation ([2])). This is because the 
Pressure model adopts a different EOS for the cold gas with 
a smaller value of 7e// = 4/3 than the usual 7 = 5/3. As we 
mentioned in Section [2.21 the value of 7e// = 4/3 has other 
benefits that it can prevent the artificial fragmentation of 
gas. Although the slope is slightly shallow, it is still in the 
acceptable range of current observational uncertainties. 



5 COSMIC STAR FORMATION HISTORY 
AND GALAXY EVOLUTION IN THE NEW 
STAR FORMATION MODEL 

In this Section, we compare the results of the SH model and 
the Pressure model on the cosmic SF history and galaxy 
evolution. Figure [5] shows the cosmic SF history for both 
models. Here, in order to cover the entire history, we in- 
clude the results from a larger volume but lower resolution 
simulation, which are initially made up of 2 x 400'^ gas and 
dark matter particles in a 100/i~^ Mpc comoving box down 
to z = (the N400L100 series). The combination of N216L10 
and N400L100 series allows us to cover both high and low 
redshifts, and alleviate any resolution effects. 

At early times, the Pressure model shows lower SFR 
than the SH model. Since the two simulations have the same 
cosmology and the same initial conditions, the dynamical 
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evolution of gas is identical until the nonlinear growth and 
star formation is initiated. As we discussed in Sections 12.21 
and |4l the star formation is suppressed in the low-density 
regions in the Pressure model compared to the SH model, 
and it also has a higher density threshold for star formation, 
resulting in the suppression of star formation at early times. 

5.1 Peak Redshift 

More interestingly, the suppression of early star formation 
shifts the peak of p-i,[z) to a lower redshift. In the Pressure 
model, the peak moves to a lower redshift by /S.z = 1 — 2 
compared to the one in the SH model. 

The total amount of stellar masses at the end of the sim- 
ulations are similar in the two models owing to the same cos- 
mology and the same structures. The SH model forms more 
stars at high-z, but the SF slows down earlier because the 
total amount of cold gas is limited. In contrast, the Pressure 
model suppresses the early star formation and leaves more 
cold gas available at lower redshifts, therefore the SFR slows 
down later than in the SH model. The shift of the peak in 
p*(z) results from the combined effect of suppression of the 
early star formation and the limited amount of cold gas at 
late times. Based on Figure (5] we expect that the true peak 
is located in-between the peaks of N216L10 and N400L100 
runs, i.e, at around 2 = 2 — 4. 

The location of the peak redshift of the cosmic SFR 
has shown co nsiderable discrepancy betwee n observations 
and theories. ISpringel fc Hernguistl (|2003bl ) claimed that 
the peak of Pi,{z) lies earlier than 2 = 5, while obser- 
vations suggest that th e peak is at 2 5-, z < 4 (e.g., 
iHopkins fcBeacomll2006l ). [Nagamine et al.l (|2004l ') compared 
two different types of hydrodynamic simulations, SPH and 
Eulerian TVD codes, and found a good agreement in f)-,,{z) 
between the two simulations with the peak being at z > 4. 
iNaeamine et al ] I2OO&I also compared the predicted stellar 
mass densities p* with observations, and found that theory 
predicts higher than the observational estimates. They 
both concluded that the early SF seems to be a generic fea- 
ture of the concordance ACDM model, and suggested that 
the current observations could be missing nearly half of the 
Pi, in the Universe at high-z. 

In Figure [21 we also compa re our simulation results with 
the theoretical fitting model of iHernguist fc Springe! l|2003l . 
the HS model), which is based on a series of cosmological 
SPH simulations combined with semi-analytic arguments. 
Since our simulations include metal cooling, the star forma- 
tion is enhanced by 20 — 50% co mpared to the simulations 
without metal cooling (jChoi fc N agamine 2009). To include 
this enhancement, we multiply a factor of 1.3 to the original 
HS model formula in Figure [S] As expected, the HS fitting 
agrees well with the SH model. 

We also show two compilations of observational es- 
timates of p ^ jz) i n Figure [S] The cyan shading is from 
iKistler et all (|2009l ). which combines the data inferred from 
the Swift gamma-ray b ursts (GRBs), the data compiled 
by iHopkins fc BeacomI (l2006f). and the high-z UV data 
from 'B ouwens et~al] (|2007l . 1200^ 1. The high-z SFRs in- 
ferred from GRBs tend to be higher than the previous esti- 
mates, which makes the cyan shading quite wide. The yel- 
low shading shows the observational estimates complied in 
iNagamine et al.l (|2006l 'l. who used different dust corrections 



from iHopk ms fc BeacomI (|2006l '). Owing to different dust 
corrections, the yellow shading is lower at \ow-z, and shows 
better agreement with our simulations. The current obser- 
vations of p*(z) still show large uncertainties particularly 
at high-2, and seem to suggest that the peak of pi,{z) lies 
at z < 4, which is at a lower redshift than the SH model 
prediction. 

The peak of Pi,{z) will be better determined in the 
near future using the data from the next generation of tele- 
scopes, such as the 30 meter telescopes and the James Webb 
Space Telescope {JWST). These telescopes will measure the 
faint-end slope of the luminosity function more accurately, 
and the estimates of the UV luminosity density from star- 
forming galaxies will become more accurate. Our current 
results suggest that the adopted SF model was partially re- 
sponsible for the earl y peak of p-,,[z) in the simulations of 
ISpringel fc Hernguistl (|2P03b) and the HS model. We find 
that the imp roved SF mo del can mitigate the conflict that 
was found bv lNagamine et al.l (|2004l ) as we describe below. 

5.2 Stellar Mass Density 

The delayed star formation in the Pressure model gives rise 
to changes in the other observables as well. Figure [6] shows 
the evolution of the global stellar mass density p* as a func- 
tion of redshift. Several observations h ave es timated the evo- 

'^'2003'; 



lution of p-k(z) (e.g., Dic kinson et al 



. iRudnick fc etal 

I2OO6I : iPerez-Gonzalez et al..>200a : .Marchesini et al.ll2008l ')Tlri 
thi s figure, we show the ran ge of observational data compiled 
bvlMar chesini et al.l (j2008l ) with a yellow shading, obtained 
by integrating the observed stellar mass function. Note that 
they fit the observed data at > 10^° with a Schechter 
function by a maximum likelihood method, and the integra- 
tion down to A/t = 10* Af© is based on the extrapolated 
Schechter mass function fit. Compared to these observa- 
tional estimates, our simulations predict somewhat higher 
p*. Since the p* in the Pressure model tends to be lower 
than in the SH model, it shows better agreement with the 
current observational estimates. 

It is still possible that the current observations miss 
a considera ble number of low-mas s galaxies at high-z, as 
discussed bv lNagamine et al] l|2004h . To quantify the missed 
fraction of p*, we plot p*(2) with and without the mass- 
cut of Mi, > 10* Af© in Figure |6]d. Comparison of the two 
estimates shows that the measurement with the mass-cut 
underestimates p* particularly at high-2, e.g., by ~65% at 
2 = 6. This difference is evident only in the N216L10 series, 
whose galaxy resolution reaches < 10^ Mq. In the N400L100 
series, the galaxy resolution is above W^Mq, therefore the 
discrepancy between the two estimates does not show up 
in Figure [SJs. Since most of current p* is computed with 
galaxies with Mi, > 10^ Mq, they could be missing as much 
as 65% at 2 = 6 and 30% at 2 = 3. If we change the mass 
limit to Mi, = lO^^Af©, which is the flux limit of current 
high-2 galaxy surveys, instead of IO^A/q, then the missed 
stellar mass fraction would be more than 50% at 2 = 3. 

Finally, we point out that there is an interesting contra- 
diction between Figure [5] and Figure [6] The observed pi,{z) 
tends to be higher than our simulations except at very high- 
2, while the observed p*(2) tend to be lower than the sim- 
ulations. This c o ntrad iction, which was also discussed by 
iNagamine et al.l l|200d) . may result from the missed low- 
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(a) All stars (b) Stars in the galaxies with A/* > 10* Mq 

Figure 6. Evolution of the global stellar mass density in the simulations as a function of redshift. Left panel: The stellar mass 
includes all the star particles in the simulations. Right panel: The stellar mass is estimated by integrating the stellar mass function only 
over the range of Af*^ > lO^AfQ. The solid lines are for the fiducial run of N216L10 series, and the dashed lines are for the N400L100 
series. The results of the two models are shown (the SH model in red and the Pressure model in blue) for each simulation. The yellow 
shading represents the range of observational estimates shown in lMarchesini et al. I bOOSl ). The SH model predicts higher p* owing to the 
early star formation. The estimate with > 10^ Mq mass-cut (panel b) misses the stellar mass density at high redshift. 



mass galaxies at high-z or the redshift evolution of the 
stellar initial mass function l|Davell2008l : Ivan Dokkumll2008l : 
IWilkins et al.! 2008). Future observations by 30 m telescopes 
and JWST may provide an answer to this conundrum. 



5.3 Baryonic Mass Functions 

Different cosmic star formation histories may result in differ- 
ent galaxy evolution. We compare the baryonic mass func- 
tions for two different SF models in Figure ^ At z — 6, 
the stellar mass functions (left column) for the two mod- 
els show noticeable differences: the galaxies in the Pressure 
model tend to have lower stellar masses (or the number of 
galaxies for a given is lower). 

However, it does not mean that the total baryonic mass 
of galaxies is lower in the Pressure model. The galaxy bary- 
onic mass functions (right column) show that the galaxies 
in the Pressure model run are sometimes even more mas- 
sive than those in the SH model run. This is because the 
star formation is suppressed in the low density regions in 
the Pressure model run with a higher SF threshold density, 
which results in lower SFRs at early times compared to the 
SH model run. The subsequent supernova feedback is also 
weaker in the Pressure model, leading to more efficient gas 
accretion. Consequently, the galaxies in the Pressure model 
run can have higher baryonic masses than those in the SH 
model at early times. 

The difference seen at z=6 in the stellar mass functions 
between the two SF models decreases at low redshift. At 
2 = 3, the two SF models show similar stellar mass functions 



except at the very massive end. The SH model has a larger 
number of massive galaxies than the Pressure model, and 
agrees better with the observational data (cyan and yellow 
shadings) . The larger number of massive galaxies in the SH 
model can be ascribed to the merger of smaller galaxies that 
formed earlier more efficiently than in the Pressure model. 
The Pressure model cuts across the observed range, and is 
at the higher and lower edge of the observed range at M* ~ 
10^°Mq and lO^'^M©, respectively. 

At z — 0, the two SF models show similar mass func- 
tions. Both models agree well with the observation at A/* < 
lO^^Af0, but they both overpredict the observed data signif- 
icantly at the massive end. This discrepancy may be due to 
the lack of AGN feedback in our simulations, which is con- 
sidered to be the major mechanism to quenc h star formation 
in massive galaxies (e.g.. lCroton et aL I I2OO6I ). We plan to m- 
clude the evolution of supermassive black holes and their 
feedback effects in our future simulations. 



5.4 Gas Fraction 

The suppression of early star formation increases the gas 
fraction in galaxies. Figure |8] shows the mean gas fraction 
of galaxies as a function of galaxy stellar mass for the two 
SF models. The mean is defined as the ratio of the total 
gas mass to the total baryon mass (gas-|-stars) for all the 
galaxies in each mass bin, i.e., J2i -^gas.i/Zli Mbaryon,i. The 
galaxies in the Pressure model run are more gas rich than 
those in the SH model run at all redshifts, except for massive 
galaxies at, z — 0. The balance between galactic outflow 
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Figure 7. Baryonic mass functions for the two star formation models: the SH model (solid lines) and the Pressure model (dashed lines). 
The left column panels show the stellar mass functions, and the right column panels show the baryonic (gas + star) mass functions. 
At z = % and 2 = 3, in each panel we plot three different sets of simulations with different resolution to cover a wide mass range: the 
N216L10 runs cover the low-mass end, the N400L34 runs cover the intermediate mass range, and the N400L100 runs cover the most 
massive end of the mass functions. At 2 = 0, we only show the N400L100 run, because the runs with small volumes are stopped at higher 
redshifts. The shade d regions at 2 = 3 repres ent the range of observed stellar mass functions at 3 < 2: < 4 (yellow shading) an d 2 < z < 3 
cyan shading) from lMarchesini et al.l 1I2OO8I ). The green shading at 2 = shows the observed local stellar mass function from lCole et al.l 
200 ih . The two SF models show different mass functions at 2 = 6, but this discrepancy reduces at 2 = 3 and 2 = 0. 



and gas accretion determines the gas fraction of the massive 
galaxies at 2 = 0. We will discuss the eflects of different 
galactic outflow models in a separate publication. Overall, 
the gas fraction decreases with decreasing redshift. 

One might expect that the increased gas fraction in 
the Pressure model may enhance the star formation asso- 
ciated with the mergers of gas-rich spiral galaxies. Recent 
simulations and models show that the SFR originating from 
starburst during gas -rich galaxy mergers is ~ 10% of the 
total spheoid mas s teox e t al.1 12OO8I : iHopkins et al.1 l2009l : 
iHopkins fc Herriguist 2010 ). Therefore, we expect that the 



increase of gas fractions in the Pressure model would not 
change the total amount of star formation from mergers in 
our simulations very much. 



6 SUMMARY 

We have used the cosmological N-body/Hydrodynamics sim- 
ulation to study the effects of different star formation models 
on the cosmic star formation history and galaxy evolution. 
Our main results are as follows; 
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Figure 8. The mean gas fraction of galaxies as a function of galaxy stellar mass. The galaxies in the Pressure model run are more gas 
rich than those in the SH model run at all redshifts, except for the most massive galaxies at z = 0. 



(i) The change of cosmological parameters alters the am- 
plitude of the cosmic SFR density, but does not change the 
shape of the cosmic SF history very much. The cosmolog- 
ical parameters related to the primordial power spectrum 
(ris and erg) affect early star formation history significantly, 
while those related to the matter contents {Q,m, ^a, and Slj) 
effect the SF history more at low redshift with a smaller 
degree. We find that the cosmological parameters hardly 
change the local SFR for a given H I column density. 

(ii) We developed two new star formation models, which 
consider the effects of H2. Both the Blitz model and the Pres- 
sure model reduce the overprediction of projected SFR at 
low A^Hi. However, the Pressure model gives a more realistic 
SF threshold density and Ssfr, for a given A^hi than the Blitz 
model, in better agreement with the empirical Kennicutt- 



Schmidt law. Therefore, we treat the Pressure model to be 
our new fiducial SF model. We will continue to refine our 
SF model in the future. 

(iii) The Pressure model predicts a large number of 
columns with log Esfr < —4.5, which are below the thresh- 
old of the current observations. If our result is correct, then 
the current surveys o f nearby spiral galaxies (e.g., SINGS 
iKennicutt et al.ll2003l ') could be missing ~0.005% of the to- 
tal star formation and ~30% of columns with these low SF 
surface densities. 

(iv) The Pressure model reduces the SFR in low-density 
regions, which causes the suppression of early star formation. 
Owing to this suppression, the peak of the cosmic SF history 
is shifted to a lower redshift, making our results to be more 
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consistent with the recent observational estimates. The shift 
of the peak also decreases the global stellar mass density. 
We find that the Pressure model still predicts higher stellar 
mass density than the current observational estimates at 
high redshift. If our results are correct, then the current 
observations could be missing as much as 65% and 35% of 
total stellar masses at z=6 and 2=3, respectively, due to 
the flux limit of the surveys. In particular, our simulations 
contain a large population of low-mass galaxies with Af* < 
10* Mq at z > 3 that are undetected by the current surveys. 

(v) We point out an interesting inconsistency between the 
observational estimates of the cosmic SFR density and the 
global stellar mass density (Figure [5] vs. Figure 16]). While 
the simulation results of the SFR density is on the lower 
side of the observed range, the predicted stellar mass density 
is larger than the observed range. This suggests that there 
is an inconsistency between the observational estimates of 
the two quantities. The only uncertainty in this argument is 
the stellar IMF, which changes the amount of recycled gas 
and stellar luminosity output per unit mass of collapsing gas 
cloud. Some researchers have invoked a top-heavy IMF to 
solve this problem (see the discussion in Section [SJ. 

(vi) Owing to the suppression of the early star formation, 
the galaxies in the Pressure model run tend to have lower 
stellar masses, but similar total baryonic masses to those in 
the SH model. Therefore, the galaxies in the Pressure model 
run tend to be more gas rich than those in the SH model. 



In this paper, we mainly focused on the results from 
the N216L10 and N400L100 series simulations. In the fu- 
ture, we will carry out a larger number of simulation series 
with different resolution and different box sizes to handle 
the resolution effects better. However, we expect that our 
findings will not change very much, because the qualitative 
changes due to cosmological parameters and the SF mod- 
els that we discussed in this paper should be robust against 
resolution changes. 

We have developed new SF models, however, our treat- 
ment still does not explicitly incorporate the formation of H2 
owing to the resolution limitation. Recently, a few groups 
attempted to carry out a very high resolution simulation 
with radiative transfer and explicit treatment of H2 for- 
mation and destru ction (e.g., iRobertson fc Kravtsovl [20081 : 
iGnedin et aLll2009l ). So far, they have been able to imple- 
ment this approach only for a single galaxy, but such a sim- 
ulation certainly provides importa nt physical insights for 
an im proved star formation model. [Robertson fc Kravtsovl 
confirmed that Esf r shows bette r corre lation with 
Eho , rather than with E 

gas ■ 

iGnedin et aD ()2009l ) found that 
the transition from atomic to fully molecular phase depends 
on the metallicity. It implies that the threshold density 
should be a function of gas metallicity in an improved star 
formation model. This causes very low star formation effi- 
ciency in the low-mass, low-metallicity galaxies, and one can 
expect further suppression of star formation in high-z galax- 
ies that are not enriched with metals yet. In the future, we 
will attempt to incorporate the effects of H2 more explicitly 
and improve our SF models by considering the metallicity 
effects in cosmological simulations. 
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